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ABSTRACT 

We conduct a search for z > 7 dropout galaxies behind 11 massive lensing clusters using 21 arcmin^ 
of deep HST NICMOS, ACS, and WFPC2 image data. In total, over this entire area, we find only 
one robust z ^ 1 z-dropout candidate (previously reported around Abell 1689). Four less robust z- 
dropout and J-dropout candidates are also found. The nature of the four weaker candidates could not 
be precisely determined due to the limited depth of the available optical data, but detailed simulations 
suggest that all four could be low-redshift interlopers. We compare these numbers with what we might 
expect using the z ^ 1 UV luminosity function (LF) determined from field searches. We predict 2.7 
z ~ 7 z-dropouts and 0.3 z ~ 9 J-dropouts over our cluster search area, in reasonable agreement with 
our observational results, given the small numbers. The number of z > 7 candidates we find in the 
present search are much lower than has been reported in several previous studies of the prevalence of 
z > 7 galaxies behind lensing clusters. To understand these differences, we examined z > 7 candidates 
in other studies and conclude that only a small fraction are likely to be z > 7 galaxies. Our findings 
support models that show that gravitational lensing from clusters is of the most value for detecting 
galaxies at magnitudes brighter than L* {H < 27) where the LF is expected to be very steep. Use of 
these clusters to constrain the faint-end slope or determine the full LF is likely of less value due to the 
shallower effective slope measured for the LF at fainter magnitudes, as well as significant uncertainties 
introduced from modelling both the gravitational lensing and incompleteness. 
Subject headings: galaxies: evolution — galaxies: high-redshift 



1. INTRODUCTION 

Because of the great distances and extreme faintness 
of galaxies at z > 7, as well as the high sky backgrounds, 
detection of galaxies at such high redshifts remains ex- 
tremely challenging. It is not surprising that the number 
of robust z>7 candidates is still very small (see, e.g., 
Bouwens & Illingworth 2006; Bouwens et al. 2008; Oesch 
et al. 2008). Gravitational lensing by galaxy clusters has 
been highlighted as an efficient way of improving this 
situation, due to the significant areas on the sky behind 
these clusters with sizeable magnification factors to am- 
plify light from faint sources. However, this advantage is 
offset by the greatly reduced source plane volume in the 
highly magnified regions. 

Because of the trade-off between depth and area, the 
utility of clusters for these searches depends strongly 
on the slope of the luminosity function (e.g., Broad- 
hurst et al. 1995). If the effective slope of the LF 
{—d{\ogd^)/d\ogL) is greater than 1, the gain in depth 
more than compensates for the loss in area, increasing 
the overall number of sources (e.g., Broadhurst et al. 
1995) over that found in the field. Such steep slopes 
are found at magnitudes brightward of L* (correspond- 
ing to H<27 AB mag for z ~ 7 galaxies). By contrast, 
at fainter magnitudes (H>27 for z ~ 7 galaxies) the ef- 
fective slope of the LF is not quite so steep (e.g., the 



faint-end slope for the Bouwens et al. 2007 z ~ 6 LF 
corresponds to -(i(log$)/d(logL) 0.7 < 1). This 
trade-off between depth and area is such that the sur- 
face density of dropouts is lower behind clusters than in 
the field at faint magnitudes. Overall these considera- 
tions suggest that the most significant advantages will 
be achieved at bright magnitudes (>L*) where the LF 
is very steep. Shallow searches over many clusters, in 
particular, would seem to be the most rewarding. 

Here we assess the promise of clusters for studying 
z > 7 galaxies by conducting a careful search for high 
redshift galaxies in all the currently available HST NIC- 
MOS imaging data over 11 massive low-redshift galaxy 
clusters. Bradley et al. (2008) have already conducted 
such a search around Abell 1689 and reported one highly 
robust z > 7 galaxy. Richard et al. (2006) examined 2 
clusters and reported 13 z > 6 candidates, while Richard 
et al. (2008) have examined 6 clusters and reported 12 
z > 6 candidates. The present paper represents an in- 
dependent assessment of the prevalence of these sources 
behind massive low-redshift clusters. We take advan- 
tage of ~21 arcmin^ of very deep NICMOS data behind 
11 lensing clusters with optical ACS-I-WFPC2 coverage, 
7 of which were already considered in the Bradley et 
al. (2008) and Richard et al. (2008) papers. For our 
z > 7 search, we will utilize many of the same photomet- 
ric techniques we have employed over the past few years 
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Fig. 1. — Optical ACS images of the 11 massive lensing clusters we used to search for star-forming galaxies at z > 7. The name of 
each cluster is included in the upper left corner of the corresponding cluster image. Overlaid on these images are outlines of our deep 
NIC3 search fields, where magenta indicates that the field has deep (~3-5 orbit) data in both the Jiio and Hieo bands, blue indicates the 
field only has shallow (~l-2 orbit) data in the Jiio band, and red indicates the field only has deep (~4 orbit) data in the -ffieo band (see 
Table 1). The green contours show the position of the critical curves at (i.e., where the magnification /i>100) as determined from 

several published lensing models (Limousin et al. 2007; Limousin et al. 2008; Elfasdottir et al. 2007; Jee et al. 2007; Franx et al. 1997). 



TABLE 1 

HST NICMOS DATA AROUND MASSIVE GALAXY CLUSTERS USED TO SEARCH 
FOR Z > 7 GALAXIES. 



Area So- Depth'' NICMOS 



Name 


(arcmin^) 


Optical'' 


Z8l 


w 


-^110 




orbits 


RoF 


MS1358 


1.4 


28.0 


27, 


,5 


26.8 


26.7 


18 


[1] 


CL0024 


1.4 


27.8 


27, 


,4 


26.8 


26.7 


18 


[1] 


Abell 2218 


1.2 


27.9 


27, 


,6 


26.8 


26.7 


17 


[1] 


Abell 2219 


1.4 


27.5 


27, 


,2 


26.8 


26.7 


18 


[1] 


Abell 2390 


1.5 


26.6'* 


26, 


,9 


26.8 


26.7 


18 


[1] 


Abell 2667 


1.5 


26.5'* 


26, 


,9 


26.8 


26.7 


18 


[1] 


Abell 1689 


5.7 


28.1 


26, 


,8 


26.4 




19 


[2] 


Abell 1703 


3.1 


27.9 


26, 


,7 


26.0 




6 




lE0657-56<= 


2.7 


27.6 


27, 


,2 


26.0 




4 




Abell 1835 


0.7 


27.1<* 


27, 


,2 




26.7 


4 




AC114 


0.7 


27.1<* 


27, 


,2 




26.7 


4 




Total 


21.3 












144 





''StT limits assume a 0.3"-diameter aperture for ACS/WFPC2 and 0.6"- 
diameter aperture for NICMOS. 

''The depth of the deepest single-band optical (< O.Sfi) image available over 
the cluster. Note that some of the clusters listed here have very deep data 
in > 2 bands (particularly MS1358, CL0024, Abell 2218, Abell 1689, Abell 
1703), so the effective depth of the combined optical data is often >0.4 mag 
deeper than tabulated here. 

'^References: [1] Richard et al. 2008, [2] Bradley et al. 2008. 

''The optical data available over these four clusters are from WFPC2 and are 
only moderately deep. As a result, we might expect a small number of low 
redshift interlopers to make it into 2^7 selections (§4) conducted over each 
of these clusters. Not surprisingly, a substantial fraction of the weaker z ^7 
candidates in our own selection, and that of Richard et al. (2008) are found 
over these four clusters. 

<=The "Bullet" cluster 
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to identify large samples of z^4-8 BViz dropouts in the 
field (Bouwens et al. 2006: Bouwens et al. 2007; Bouwens 
et al. 2004; Bouwens et al. 2008). 

We begin this paper with a summary of the observa- 
tional data available to search for z > 7 galaxies behind 
galaxy clusters (§2). In §3, we describe our techniques for 
constructing source catalogs and doing the photometry 
(§3). In §4, we summarize the results of our selection and 
discuss sources of contamination. In §5, wc compare our 
results with (1) what we would expect based upon the 
2; ~ 7 LF derived in the field and (2) the search results 
from other teams. Wc then conclude by discussing the 
implications of this study (§6). In the Appendix, we pro- 
vide an assessment of the z > 7 candidates reported by 
Richard et al. (2008). In particular, we assess the mer- 
its of searching for z > 7 galaxies behind clusters versus 
searching for these galaxies in the field (§6). Wc assume 
r^o = 0.3, Oa = 0.7, Ho ^ 70km/s/Mpc throughout. Al- 
though these parameters are slightly different from those 
determined from the WMAP five-year results (Dunkley 
et al. 2008), they allow for convenient comparison with 
other recent results expressed in a similar manner. The 
HST filters F555W, F606W, F625W, F702W, F775W, 
F814W, F850LP, FllOW, and F160W will be denoted as 

^555, ^06, ''6255 -^702, 1775, hli, ^850, ■^110, and Hiqq, 

respectively. We will express all magnitudes in the AB 
system (Oke & Gunn 1983). 

2. OBSERVATIONAL DATA 

We conduct searches for dropout galaxies over the 11 
low-redshift galaxy clusters with deep near-IR NICMOS 
and optical ACS coverage (Table 1). The near-IR data 
here were taken with the goal of finding z > 7 galaxies. 
The NICMOS coverage of the first six clusters considered 
here (CL0024, MS1358, Abell 2218, Abell 2219, Abell 
2390, Abell 2667) extends over 8.4 arcmin^ (12 NIC3 
pointings) and includes very deep imaging in both the 
Jiio and Hiea bands (Richard et al. 2008). The NIC- 
MOS coverage of three other clusters considered here 
(i.e., Abell 1689, Abell 1703, and 1E0657-56) is some- 
what shallower in general, extends over ~10 arcmin^ (17 
NIC3 pointings), and is mainly in the Jno band. The fi- 
nal two clusters considered here (Abell 1835 and AC114) 
only have deep near-IR NICMOS data in the ifieo band 
and over one ~0.8 arcmin^ NIC3 pointing per cluster. 
The layout of our NICMOS search fields is illustrated 
graphically in Figure 1. 

Deep optical imaging data are necessary for the se- 
lection of ^ > 7 galaxies and can be a significant limita- 
tion, if not available. MS1358, CL0024, Abell 1689, Abell 
1703, and Abell 2218 all possess useful optical data, with 
> 2 orbits in each of the ACS ,9475, ?'625, and 1775 bands 
and > 6 orbits in the ACS zsso band. The optical cov- 
erage of Abell 2219, Abell 2390, Abell 2667, 1E0657-56, 
Abell 1835, and AC114 is generally shallower in depth 
and primarily with WFPC2, though ~4 orbits of ACS 
2;-band coverage are available for each. 

All the available ACS and NICMOS data over these 
clusters arc; processed into image mosaics using the ACS 
GTO pipeline "apsis" (Blakeslee et al. 2003) and NIC- 
MOS pipeline "nicred.py" (Magee et al. 2007). Reduc- 
tions of the WFPC2 data are obtained from the Cana- 
dian Astronomy Data Centre. All reductions are regis- 
tered onto the same frame as the NICMOS data. 



3. SAMPLE CONSTRUCTION 

(a) Catalog Generation: Our procedure for generating 

source catalogs and doing photometry is identical to that 
performed in Bouwens et al. (2008). Briefly, we run SEx- 
tractor (Bertin & Arnouts 1996) in double-image mode 
to do object detection and photometry. For the detection 
image, we use the square root of the image (Szalay 
et al. 1999), which we construct from the NICMOS Jno 
and HiQo images for our zs5o dropout selection and the 
NICMOS i?i6o band image for our Jno dropout selec- 
tion. Our photometry is then conducted on our ACS, 
WFPC2, and NICMOS images, which are point-spread 
function matched to the NICMOS -ffieo image. Colors 
are measured in small-scalable (Kron 1980) apertures, 
assuming a Kron (1980) factor of 1.2. These fluxes are 
then corrected to total magnitudes using the light within 
a larger Kron (1980) aperture (adopting a Kron factor of 
2.5). These latter corrections are made using the square 
root of the image to improve the S/N. Figure 5 of 
Coe et al. (2006) provides a graphical description of a 
similar multi-stage procedure for measuring colors and 
total magnitudes. Typical aperture diameters are 0.3" 
and 0.6" for color and total magnitude measurements, 
respectively. 

(b) Selection Criteria: We use the same selection 
criteria for identifying star-forming galaxies at 2; > 7 
that we used in our previous work on the identiflca- 
tion of such galaxies in field data sets like GOODS or 
the HUDF (e.g., Bouwens et al. 2008). Speciflcally, 
WC! rcqiiire a z-dropout candidate to satisfy the crite- 
rion ((z850 - Jiia)AB) > 0.8) A ((z85o - Jiw)ab > 
0.8 + 0.4( Jno — -ffi6o)AB) where A represents the logi- 
cal AND symbol. J-dropout candidates are expected 
to satisfy the criterion (Jno — Hiqo)ab > 1-3. In cases 
of a non-detection in the dropout band, we set the flux 
in the dropout band to be equal to the Icr upper limit. 
We require each candidate to be completely undetected 
(< 2a) in all passbands blueward of the dropout band. 
We also demand that each candidate be detected at 5.5cr 
in the iJieo band in a 0.5"-diameter aperture to eliminate 
spurious sources. 

Our selection criteria are modified slightly for clusters 
that do not have NICMOS imaging in both the J and H 
bands. For clusters with only J band imaging, we apply a 
(-2^850 ~ Jiio)ab > 1.0 criterion to select z^7 z-dropouts 
and for clusters with only H band imaging, we apply 
a (^850 — Hiqo)ab > 1-2 criterion to identify possible 
star-forming galaxies at z > 7. Both criteria should be 
successful in identifying candidate z>7 galaxies (albeit 
with a higher contamination level than selections relying 
on both J and H data). 

4. RESULTS 

After careful application of our selection criteria to all 
11 clusters under study here, we identify 4 z-dropout 
candidates and 1 J-dropout candidate. We also un- 
cover a small number of candidates that appeared to 
be promising z > 7 candidates (e.g., the candidates at 
00:26:37.90, 17:09:10.4 or 00:26:35.11, 17:10:10.3 behind 
CL0024), but which show modest (~2a) detections in 
passbands blueward of ^850 and therefore were excluded. 

Postage stamps of the four candidate star-forming 
galaxies at z > 7 are provided in Figure 2. Other prop- 
erties are given in Table 2. Our candidate in Abell 1689 
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TABLE 2 

Z>7z, J-DROPOUT CANDIDATES. 



Object ID 


R.A. 


Dec 


0.6/i - Jiio 


^ 0.8m - Jiio'' 


2850 — J^llO 


JllO — Hl60 












7 z-dropouts 








Aiesg-zDit: 


13:11:29.73 


-01:19:20.9 


>2.5'l 


>2.5<^ 


>2.2d 


0.6±0.2 


24.6±0.1 


A2390-zDl'= 


21:53:34.09 


17:41:41.1 


>1.7'^ 


>1.5<l 


1.1±0.8 


0.8±0.2 


25.2±0.2 


A2667-zDl'= 


23:51:40.06 


-26:05:13.8 


>0.9'' 


>0.4d 


>1.2<^ 


0.0±0.3 


26.1±0.2 


A2667-zD2'= 


23:51:36.85 


-26:05:21.4 






0.9±0.4 


0.0±0.3 


25.6±0.2 








z ~ 


9 J-dropouts 








A2390-JD1'='S 


21:53:34.12 


17:41:44.2 


>1.2d.f 


>l.ld.t 


>1.3'i'' 


>1.9 


26.0±0.2f5 



For uniformity of our analysis, all of the 2 > 7 candidates here are from our own z > 7 z and J dropout searches in 
cluster fields and do not include independent search results. 
^This colour corresponds to r625 — Juo for Abell 1689, V555 — Jiio for Abell 2390, and Veoe — ^110 for Abell 2667. 
''This colour corresponds to 1775 — Jno for Abell 1689 and 7814 — •.'iio for Abell 2390 and Abell 2667. 
'^A1689-zDl is our most robust 2^7 candidate and was previously presented in Bradley et al. (2008). 
''Lower limits on the measured colors are Icr limits. 

"No deep optical data are available for the four candidates: A2390-zDl, A2667-zDl, A2667-zD2, and A2390-JD1 (see 
Table 1). We consider all four to be relatively weak 2 > 7 2-dropout candidates. 

^The colours tabulated here are relative to the if-band, not the Jno-band, and hence are 0.6/i — H\qq, 0.8/U — Hiqq, 
and 2850 - Tfieo- 

^While this source is formally a 5.5(T detection in our reductions (and using our photometric procedures), there is some 
chance it may still be spurious. 



yI\ Ml z J H 




A2390-.TD1 



Fig. 2.— r625,«775,2850,Jii(),J7l60 images (4" X 4") of 4 2~7 
2-dropout candidates and one 2~9 J-dropout candidate (A2390- 
JDl), identified behind 11 massive lensing clusters. We consider 
the top source (A1689-zDl: Bradley et al. 2008) to be the only 
highly robust 2 > 7 candidate in our cluster search fields. The 
leftmost two images correspond to the WFPC2 Vsss and /814 bands 
for Abell 2390 and the WFPC2 V6O6 and /sm bands for Abell 2667. 
Deep WFPC2 data are not available for A2667-zD2. 

has already been reported upon before by Bradley et al. 

(2008) and is by far the strongest candidate. The SED 
derived from the very deep optical, NICMOS, and IRAC 
data can only be successfully fit by the SED of a star- 
forming galaxy at z ~ 7.6 (Bradley et al. 2008). 

Possible Contamination of our Selection by Lower- 
Redshift Galaxies: The four new sources in our selec- 
tion are much less robust z > 7 candidates. All are in 
clusters for which the optical data is much shallower. 



The 1^555, Veoe, or WFPC2 coverage over Abell 2390 
and Abell 2667 reaches to only ~ 26.6 and ~ 26.5 {5a, 
0.3" -diameter aperture) over Abell 2390 and 2667, re- 
spectively. These WFPC2 data are > 1 mag shallower 
than the ACS data available over 5 other clusters stud- 
ied here (Table 1). This will result in a much higher 
contamination rate from low-redshift galaxies scattering 
into our color selection due to noise. 

To estimate the expected contamination level, wc start 
with a selection of intermediate magnitude galaxies in the 
HUDF NICMOS field (Thompson et al. 2005) and add 
noise to the fluxes to match the errors for sources in our 
search fields. On average, we find ~ 1 low-redshift con- 
taminant over Abell 2390 and ^ 2 such sources over Abell 
2667 in these simulations. This is comparable to the ob- 
served numbers and suggests that A2390-zDl, A2667- 
zDl, A2667-zD2, and A2390-JD1 may be low-redshift 
interlopers. As an alternate estimate of the contamina- 
tion level, we start with the catalog of sources in the 
three clusters with the deepest ACS and NICMOS data 
(Abell 2218, MS1358, CL0024) and add noise (again to 
match the flux errors for sources behind Abell 2390 and 
Abell 2667). On average, we predict ~ 2 low-redshift 
contaminants over Abell 2390 and ^3.5 low-redshift con- 
taminations over Abell 2667. Again these numbers are 
comparable to the observed number of candidates over 
these clusters and suggest that all four candidates be- 
hind these clusters may be interlopers. The modest dif- 
ferences in the predicted number of contaminants for the 
two methods is consistent with what one would expect 
from small number statistics. 

Surface Density of Robust z > 7 Candidates Behind 
Massive Galaxy Clusters: We tabulate the surface den- 
sity of robust z > 7 candidates found behind clusters 
as a function of magnitudes in Table 3. Also included 
in this table are the surface densities found in the field 
(Bouwens et al. 2008; R.J. Bouwens et al. 2008). The sur- 
face density of 2 > 7 candidates behind clusters appears 
to be substantially larger than in the field at bright mag- 
nitudes (i.e., H < 25.5 or >5L*). This is exactly what we 
expect as a result of the slope of the LF at bright magni- 
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tudes. Interestingly, current observations do not provide 
any evidence for an enhancement in the surface density 
of 2 > 7 sources between H ~ 25.5 and H ~ 27.0. This is 
despite the expected steep slope of the LF at such mag- 
nitudes. This could reflect the small numbers of sources 
involved here and point to the need for imaging more 
clusters to improve the overall statistics (see also §6 and 
R.J. Bouwcns et al. 2008, in prep). 

At > 27 AB mag we are probing faintward of L* and 
likely reaching the regime where the effective slope of the 
LF is only moderately steep (i.e., — (i(log$)/(i(logL) ~ 
0.7 < 1: Bouwens ct al. 2007). For such slopes, the 
trade-off between depth and area is such that the surface 
density of z > 7 dropouts behind clusters will be lower 
than observed in the field. 

5. DISCUSSION 

5.1. Comparison with Model Expectations 

We can compare the number of > 7 galaxy candi- 
dates found in our search with that expected based on de- 
terminations of the UV luminosity function (LF) at z ~ 7 
from the field (Bouwens et al. 2008). For this calculation, 
we need Icnsing models for the clusters under study. We 
have such Icnsing models for 5 of the search clusters un- 
der study here, namely Abell 1689, Abell 1703, Abell 
2218, CL0024, MS1358. These models are described in 
Limousin ct al. (2007), Limousin et al. (2008), Elfasdottir 
et al. (2007), Jee et al. (2007), and Franx et al. (1997), 
respectively. We will scale our expectations from those 
clusters to our entire sample. 

Predicted Number of z >7 Candidates: Starting with 
the Bouwens et al. (2008) model for the sizes, UV col- 
ors, and UV LF for galaxies at z ^ 7 (Appendix B from 
that work), we generate very high-resolution mock im- 
ages (pixel size 0.02") of galaxies in a 2: 7 source plane 
and then remap these images to the image plane using 
available Icnsing models for these clusters. Wc then add 
these simulated fields to the real data and attempt to 
recover these sources using our cataloguing and selection 
procedures (§3). Generating these mock fields 20 times 
for each cluster and using our selection procedure to iden- 
tify z and J dropouts, we estimate that we would expect 
to find 2.7 z dropouts over all of our search fields. If we 
repeat this procedure for z ~ 8—10 J-dropouts assuming 
an extrapolation of the Bouwens et al. (2008) LF results 
to z ~ 9, 0.3 J dropouts are expected. We find only mod- 
est variations in the number of dropouts expected behind 
clusters with similar depths and survey areas (e.g., ~0.2 
predicted behind Abell 2218 vs. ^0.3 predicted behind 
MS1358). Similarly, use of slightly different lensing mod- 
els for the clusters only appears to have a modest effect 
on the numbers (e.g., use of an updated version of the 
Broadhurst et al. 2005 model only increases the predicted 
numbers behind Abell 1689 by ~20% over that predicted 
by the Limousin et al. 2007 model). 

Comparisons against the Number of Observed z > 7 
Candidates: The above predictions are quite consistent 
with the 1 robust z-dropout candidate, 3 other possible z- 
dropout candidates, and 1 possible J-dropoiit candidate 
discussed in §4, particularly given the small numbers and 
uncertain contamination levels. Because of the difficulty 
in interpreting the search results over the clusters with 
the shallower optical data (given the likely significant 
contamination levels), it also makes sense to restrict our 



TABLE 3 

Surface density of strong z '>7 candidates in the deep 

ne,\r-ir data neinxd cu'sters and ik the field. 







Surface 


Search Area 




Magnitude 


Density 


per candidate 


Location 


Range'' 


(arcmin"-^) 


(arcmin^)'' 



Cluster 


Hieo < 25.5 


U.UO_Q Q4 


91+101 

■"--15 


Cluster 


25.5 < Hieo < 26.5 


<0.1ld 


>9 


Cluster 


26.5 < Hieo < 27.5 


e 




Field 


H160 < 25.5 


<0.03* 


>36 


Field 


25.5 < Hieo < 26.5 


'^•'^'-0.04 




Field 


26.5 < H16O < 27.5 




2t\ 



''Without any magnification from gravitational lensing and adopt- 
ing the Bouwens et al. (2008) determination for L* at z ~ 7, 

-ffieo < 25.5 corresponds to luminosities > 5L*_7, 25.5 < -ffieo < 
26.5 coiicsponds to luminosities > 2L*_^ and < 5Z/*_.^, and 
26.5 < < 27.5 corresponds to luminosities > 0.8L*_7 and 

''This column is the reciprocal of the previous column. 

'^Bascd on the NICMOS data over all 11 galaxy cluster fields con- 
sidered in this paper 

dRascd on the NICMOS data over Abell 1835, Abell 2218, Abell 
2219, Abell 2390, Abell 2667, MS1358, CL0024, and AC114 

•^Ncar-IR data over cluster fields arc not deep enough to probe 
this regime well. Nonetheless, wc expect the surface density of 
z > 7 dropouts behind clusters to be lower than observed in the 
field, since at H > 27 AB mag we are probing faintward of L* and 
likely reaching the regime where the effective slope of the LF is only 
moderately steep (i.e., — a!(log <I>)/c!(log L) ~ 0.7 < 1: Bouwens et 
al. 2007). For such a slope, the trade-off between depth and area is 
such that the surface density of dropouts is lower behind clusters. 

'Based on the Bouwens et al. (2008) and R.J. Bouwens et al. 
(2008, in prep) search results for 2 > 7 dropouts in the field. The 
R.J. Bouwens et al. (2008, in prep) search takes advantage of more 
than 20 arcmin^ additional search area with NICMOS not consid- 
ered by Bouwens et al. (2008). 

SBased on data from the HUDF Thompson field (Thompson et al. 
2005), HUDF-NICPARl, and HUDF-NICPAR2 fields (Bouwens et 
al. 2008) 



comparison to only those clusters with the deepest data 
(i.e., excluding Abell 1835, Abell 2390, Abell 2667, and 
AC 114 from this comparison). Using only clusters with 
the deepest optical data, we predict 2.0 ^ ~ 7 z-dropouts 
and 0.2 J-dropouts in total. Again this is in reason- 
able agreement with the 1 robust z-dropout we find over 
this more restricted search area, given the small number 
statistics. For both this comparison and the previous 
one, our search results are clearly consistent with our 
predictions and therefore with the Bouwens et al. (2008) 
z ^ 7 LF derived in the field. From this exercise, it is 
also quite clear that z > 7 galaxy searches (using a small 
number of clusters) provide only a very weak constraint 
on the volume density of lower luminosity galaxies. For 
our search the uncertainty in the volume density of these 
sources is ~ 0.6 dex [factor of four] based on our sample 
size of one. 

5.2. Comparison with Previous Results 

Richard et al. (2006) reported finding 13 z > 6 candi- 
dates (lst+2nd category) in ^12 arcmin^ of ISAAC data 
behind 2 lensing clusters, while Richard et al. (2008) re- 
ported 12 other z > 7 candidates in ~9 arcmin'^ of NIC- 
MOS data behind six other clusters. 

The Richard et al. (2006) search results are signifi- 
cantly different from our results (13 z > 6 candidates 
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TABLE 4 

Number of strong z>7 candidates identified vs. the number 

of iist xicmos ()rp,it'< used for the search fsee sfi"). 



Location 




NICMOS 
Orbits 


Candidates 
per Orbit 


Orbits per 
Candidate*^ 


Cluster (Shallow'=)° 
Cluster (Deepd)f 
Field (Shallow'^)s 
Field (Deep'^)8 
Field (Both'') 


1 

3 
6 
9 


144 

~450 
~600 
~1050 


007+"-"l'^ 

U.UU(_g yQg 

007+° ''°'^ 
0.009lH°o^ 


1441S? 
1 1 7+56 



''Number of strong 2^7 candidates 

''Reciprocal of the previous column 
Search results for shallow near-IR data with 5o" depths less than 
'^27 AB mag (i.e., < 5 orbits). 27 AB mag corresponds to ~ L* at 
z ~ 7 (Bouwens et al. 2008). 

''Search results for deep near-IR data, with ha depths greater than 
~27 AB mag (i.e., > 6 orbits). 27 AB mag corresponds to ~ L* at 
2 ~ 7 (Bouwens et al. 2008). 

''AH 11 clusters considered here. 

^No very deep (i.e., reaching >27 AB mag at 5cr) near-IR data are 
available over clusters, but it is expected that the search efficiency will 
decrease significantly at fainter magnitudes since we will be probing 
faintwajrd of L* where the faint-end slope is likely only moderately 
steep (§6). 

^Based upon the NICMOS data considered in the Bouwens et al. 
(2008) and R.J. Bouwens et al. (2008, in prep) z ^ 7 searches 
''Search results for both shallow and deep near-IR data. 



to i?i60,AB ~ 25 over 12 arcmin^ vs. the one robust can- 
didate we find over 21 arcmin^). As we discuss in detail 
in Appendix C of Bouwens ct al. (2008), > 90% of their 
z > 6 candidates appear to be spurious, since none of 
the eleven z > Q candidates in their selection with sub- 
stantially deeper (-1-2 mag) NICMOS+IRAC coverage 
are significantly (>2a') detected. 

The prevalence of z>7 galaxies implied by the Richard 
et al. (2008) search results is also much greater than 
what we find in our searches. Since these candidates are 
reported over a subset of the clusters used in the cur- 
rent search, the differences are puzzling. To understand 
the possible differences, wc examined the Richard et al. 
(2008) sources in our own reductions of the same HST 
data using our photometric procedures. We find that 2 of 
their 12 candidates appear to be plausible z >7 galax- 
ies given our photometry. The other sources are likely 
contaminants (a detailed discussion of all 12 Richard et 
al. 2008 candidates can be found in Appendix A). The 
two good candidates also satisfy our z-dropout selection 
criteria (but are blended with foreground galaxies in our 
search). One of these two z > 7 candidates does not have 
particularly deep optical coverage, so there is a reason- 
able chance that it is at lower redshift. 

A comparison of our z>l sample with the Richard et 
al. (2008) z>l sample suggests that our selection may 
suffer from some incompleteness. The observed level of 
incompleteness is not surprising and is consistent with 
what we expect from blending with foreground sources 
('-~^35%). Moreover, since this same incompleteness is 
implicitly included in the predictions we make in §5.1 to 
compare with the observations, the conclusions that we 
draw based upon those comparisons should be fair.^ 

^ Even if we include the two best 2 7 candidates from the 



6. IMPLICATIONS 

The purpose of the present study is to increase the sam- 
ple of candidate z— 7 galaxies and to assess the potential 
of gravitational Icnsing by galaxy clusters to identify and 
quantify the properties of galaxies at z > 7. Lcnsing will 
increase the depth of the survey by the magnification 
factors, but decrease the search area by the same fac- 
tor. If the effective slope of the LF (— d(logrf$)/dlogL) 
is greater than 1, the gain in depth more than compen- 
sates for the loss in area, increasing the overall number of 
sources (e.g., Broadhurst et al. 1995). We would expect 
this effect to increase the numbers at bright magnitudes 
where the effective slope of the LF is very steep due to an 
apparent cut-off at the bright end (i.e., H < 27), but to 
decrease the numbers at fainter magnitudes {H > 27) 
where the faint-end slope is shallower than this (i.e., 
-(i(log(i$)/dlogL ~ 0.7 <1: e.g., Bouwens et al. 2007). 

Our simulations (§5.1) show the expected gains at the 
bright end of the LF for searches behind clusters (see Ta- 
ble 4). We expect three z > 7 galaxies over the present 
set of cluster data (144 NICMOS orbits) using the ob- 
served LF of Bouwens et al. (2008). The 48 NICMOS 
orbits/galaxy in the clusters for bright sources contrasts 
with the ~120 NICMOS orbits/galaxy needed in the field 
(nine z > 7 galaxies arc found in —1050 orbits of NIC- 
MOS data over the GOODS fields: R.J. Bouwens et al. 
2008, in prep). Nearly identical search procedures and 
selection criteria are used in both the cluster and field 
searches. Interestingly, the expected gains at bright mag- 
nitudes in clusters arc not reflected in the observational 
results above, if we only include our most robust candi- 
date (for -144 NICMOS orbits/galaxy). This could eas- 
ily arise because of small number statistics. Any gains in 
using clusters are likely to disappear at lower luminosi- 
ties faintward of the LF knee, since the slope is not steep 
enough, as noted above. 

Of course, lensing clusters can be used to poten- 
tially detect objects much fainter than in a field sam- 
ple, thereby possibly extending the LF function to fainter 
limits. Unfortunately, both small number statistics and 
the challenges of modelling clusters make this very diffi- 
cult, and likely not to be a very practical approach. This 
is because a substantial sample of objects is needed to 
faint limits to accurately determine the LF, as well as 
extremely accurate models of both lensing by the fore- 
ground cluster and incompleteness suffered by the lensed 
high-redshift population. Determining either of these lat- 
ter two quantities well (and without any systematics) is 
a great challenge. As a result, it can be difficult to even 
measure quantities like the faint-end slope of the UV LF 
at z— 4, z— 5, and z— 6 from current samples of g, r, and i 
dropouts behind lensing clusters (where the samples are 
much larger: see R.J. Bouwens et al. 2008, in prep). 

Our findings also underline the importance of having 
very deep optical data below the Lyman break for iden- 
tifying high-z dropout galaxies. Without such data, it 
is essentially impossible to distinguish bona-fide z > 7 

Richard ct al. (2008) selection in our sample, the total number of 
strong 2 > 7 candidates is still very consistent with the predictions 
we make in §5.1 from the Bouwens et al. (2008) field LFs. More- 
over, for the clusters with the deepest optical data (i.e.. excluding 
AC114, Abell 1835, AbcU 2390 and Abell 2667), only 2.2 2 > 7 
candidates are predicted (§5.1) which is quite consistent with the 
2 found (i.e., A1689-zDl and A2219-zl). 
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galaxies from the large number of low-rcdshift galaxy in- 
terlopers that may scatter into z^^l samples as a result 
of noise - as we found for the four weaker 2 > 7 dropout 
candidates in our selection (§4). Unfortunately, obtain- 
ing sufficiently deep optical data can be expensive and 
often requires > 2 — 3 x as much time as is spent obtain- 
ing the near-IR data. 

We thank Richard Ellis, Marusa Bradac. Johan 
Richard, Daniel Schaercr, and Dan Stark for useful dis- 
cussions. We thank Johan Richard for helping us to 



more precisely locate some of his team's lower S/N 
candidates in our own reductions of the available NIC- 
MOS data. We thank our referee for important feed- 
back that was valuable for improving the overall clar- 
ity of our manuscript. This research used the facilities 
of the Canadian Astronomy Data Centre operated by 
the National Research Council of Canada with the sup- 
port of the Canadian Space Agency. We acknowledge 
support from NASA grants HST-GO09803.05-A, HST- 
GO10874.04-A and NAG5-7697. 



REFERENCES 



Bartelmann, M., Steinmetz, M., & Weiss, A. 1995, A&A, 297, 1 

Bertin, E. and Arnouts, S. 1996, A&AS, 117, 39 

BlaJjeslee, J. P., Anderson, K. R., Meurer, G. R., Benitez, N., 

& Magee, D. 2003a, ASP Conf. Ser. 295: Astronomical Data 

Analysis Software and Systems XII, 12, 257 
Bouwens, R. J., et al. 2004, ApJ, 616, L79 
Bouwens, R. J., & lUingworth, G. D. 2006, Nature, 443, 189 
Bouwens, R.J., lUingworth, G.D., BlaJjeslee, J. P., & Pranx, M. 

2006, ApJ, 653, 53 
Bouwens, R. J., lUingworth, G. D., Pranx, M., & Pord, H. 2007, 

ApJ, 670, 928 

Bouwens, R.J., lUingworth, G.D., Pranx, M., Pord, H. 2008, ApJ, 

in press, arXiv:0803.0548 
Bradley, L.D., ct al. 2008, ApJ, 678, 647 

Broadhurst, T. J., Taylor, A. N., & Peacock, J. A. 1995, ApJ, 438, 

49 

Broadhurst, T., et al. 2005, ApJ, 621, 53 

Coc, D., Bem'tez, N., Sanchez, S. P., Jee, M., Bouwens, R., & Pord, 

H. 2006, AJ, 132, 926 
Dunkley, J., et al. 2008, ApJS, submitted, arXiv:0803.0586 
Eh'asdottir, A., et al. 2007, arXiv:0710.5636 

Pranx, M., lUingworth, G. D., Kelson, D. D., van Dokkum, P. G., 
& Tran, K.-V. 1997, ApJ, 486, L75 



Jee, M. J., et al. 2007, ApJ, 661, 728 

Kron, R. G. 1980, ApJS, 43, 305 

Lehnert, M. D. & Bremer, M. 2003, ApJ, 593, 630 

Limousin, M., et al. 2007, ApJ, 668, 643 

Limousin, M., et al. 2008, arXiv:0802.4292 

Magee, D. K., Bouwens, R. J., & lUingworth, G. D. 2007, 
Astronomical Data Analysis Software and Systems XVI, 376, 
261 

Meneghetti, M., Bartelmann, M., Jenkins, A., & Frenk, C. 2007, 

MNRAS, 381, 17 
Oesch, P. A., et al. 2008, ApJ, submitted, arXiv:0804.4874 
Oke, J. B., & Gunn, J. E. 1983, ApJ, 266, 713 

Richard, J., Pello, R., Schaerer, D., Le Borgne, J. -P., & Kneib, 

J.-P. 2006, A&A, 456, 861 
Richard, J., Stark, D. P., Ellis, R. S., George, M. R., Egami, 

E., Kneib, J.-R, & Smith, G. P. 2008, ApJ, submitted, 

arXiv:0803.4391 

Stanway, E. R., McMahon, R. G., & Bunker, A. J. 2005, MNRAS, 

359, 1184 

Szalay A. S., Connolly A. J., & Szokoly G. P. 1999, AJ, 117, 68 
Thompson, R. I., et al. 2005, AJ, 130, 1 



APPENDIX 

A. INDEPENDENT ASSESSMENT OF THE Z>7 CANDIDATES IN RICHARD ET AL. (2008) 

Richard et al. (2008) also conducted a search for z>7 galaxy candidates behind 6 of the 11 galaxy clusters considered 
in this study (i.e., CL0024, CL1358, Abell 2218, Abell 2219, Abell 2390, AbeU 2667). They reported finding 10 
promising 2; ~ 7 2;-dropout candidates and 2 promising 2; ~ 9 J-dropout candidates. Surprisingly, none of the 12 
candidates in the Richard et al. (2008) z > 7 sample make it into our own 2>7 dropout selections (§4) - nor do our 
2 > 7 candidates appear in their selection - even though our color selection criteria are essentially identical (and their 
search data a subset of ours). 

To understand the reason for these differences, we performed photometry on the sources using the same 0.3"- 
diameter (for ACS-I-WFPC2) and 0.6"-diameter (for NICMOS) apertures employed by Richard et al. (2008) for their 
photometry. We then corrected the magnitudes we measured from the ACS-I-WFPC2 or NICMOS data by 0.3 mag 
and 0.6 mag, respectively (again following the methodology of Richard et al. 2008). For the four Richard et al. (2008) 
candidates where light from the central cluster galaxy would contaminate the photometry (i.e., A2218-z2, A2390-z2, 
A2667-z2, A2667-jl), we explicitly fit the isophotes to the central galaxy and subtracted them off. A summary of our 
photometry of the 12 Richard et al. (2008) candidates is presented in Table Al. We present images of the 12 z>Jl 
candidates from the Richard et al. (2008) selection in Figure 1. 

Of the 12 candidates, two (i.e., A2219-zl and A2667-z2) seem like plausible z > 7 galaxies (A2219-zl being the 
stronger of the two). Both candidates satisfy our selection criteria (but were blended with foreground galaxies in our 
catalogs and therefore not included in our candidate lists). The other ten candidates appear unlikely to correspond to 
high-redshift sources, given our photometry (although in the case of A2667-zl it is difficult to rule out the source being 
at 2; > 7). Ascertaining the precise redshift of A2667-z2 is diSicult given the limited depth of the optical data around 
Abell 2667. Without deeper optical data, it is essentially impossible to know which 2; ^ 7 candidates correspond 
to high-redshift galaxies and which arc lower redshift contaminants. Note that we encountered similar difficulties in 
ascertaining the nature of the 4 weaker z >7 candidates foiuid in our selection behind Abell 2390 and Abell 2667 
(neither of which has particularly deep optical coverage). Our simulations and tests (§4) suggested that substantial 
contamination by low redshift galaxies is quite likely. 

In addition to the above comments about specific 2: > 7 candidates in the Richard et al. (2008) selection, we also have 
several general concerns about the properties of this zsso-dropout selection. Our first concern regards the Jno — i?i60 
color distribution of this selection. The colors are ~ 0.5 mag bluer on average than what we find in the field for our 
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Fig. 1. — r(i25,i775,Z850,'^li0,^^l60 images (4" x 4") of the 12 2:>7 candidates in the Richard et al. (2008) selection. The images are 
scaled for display in the same way as we used in Figure 2 to present the z ^7 candidates we identified behind clusters (and as used in 
presenting z > 7 candidates in the HUDF/GOODS fields: Figure 3 from Bouwens et al. 2008). The leftmost two images correspond to 
the WFPC2 Vb5B and 7814 bands for Abell 2390, the WFPC2 Rro2 band for Abell 2219, and the WFPC2 Veoe and 7814 bands for Abell 
2667. Photometry for these candidates is presented in Table Al. A2219-zl appears to be a blend of two sources, the upper right one of 
which appears to be a potential z >7 candidate. A2667-z2 satisfies our z-dropout selection criteria, but it is difficult to determine if it is 
a probable z ^7 galaxy given the limited depth of the optical data (Table 1: see also §4). A2390-zl and A2390-z2 have z — J colours that 
appear to be too blue (see Table Al) for them to be x ~ 7 z-dropouts. The colors we measure for A2667-zl do not satisfy our z-dropout 
selection criteria, but we cannot completely rule out this source as a z > 7 candidate. CL1358-zl, CL1358-z2, and CL1358-z3 are not 
detected at > Scr significance in our reductions of the NICMOS data (the possible sources seen in the J-band data have a significance of 
only 1 — 2.7(7 in our reductions). A2218-zl is sufficiently close to a bright source, as to make unambiguous detection difficult. A2218-z2 
is only found after subtracting the wings of the central cluster galaxy and is detected at only modest significance (~ 2 — 3it); the nature 
of this source is clearly very uncertain. Both A2219-jl and A2667-jl show significant detections at optical wavelengths and are almost 
certainly low rcdshift contaminants. 

z-dropout sample (Bouwens et al. 2008: Figure Al). They are also ^ 0.5 mag bluer than what one would expect for a 
young (3 = —2 star-forming galaxy population (as observed at z 5 — 6: e.g., Lchnert & Bremer 2003; Stanway et al. 
2005; Bouwens et al. 2006). It is hard to understand why z-dropouts from the Richard et al. (2008) selection would 
be so much bluer than these expectations unless the candidates had strong Lya emission in the Jno band. However, 
such an explanation would appear to be ruled out by the follow-up spectroscopy that Richard et al. (2008) conduct 
that find no such emission for the 7 z-dropout candidates they observe. 

Second, we find it worrisome that almost all (8 out of 10) of their zsso-dropout candidates only satisfy their zsso — ^iio 
color criterion by a small margin (< 0.3 mag) whereas in the Bouwens et al. (2008) field sample, almost half of the 
z-dropouts satisfy the Zs5o — Jno color criterion by a wide margin (i.e., > 1 mag). One would expect the situation 
to be quite opposite here, as a result of the much flatter number counts expected behind galaxy clusters (the number 
of bright sources in cluster fields would increase due to magnification by the cluster and the number of faint sources 
would decrease because of the loss of area) . This should result in a larger fraction of sources satisfying the selection 
criterion by a wide margin, not a smaller fraction. 



TABLE Al 

Independent photometry on the Richard et al. (2008) z > 7 candidates. 



0.6"-diameter aperture'^ 



^^DjeCt LLJ 


U.D/i— Jiio 


U.»/i— Jiio 


2860 — -.^110 




-H16O 
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' 7 2-dropouts 








CL1358-Z1'* 


>1.9 


>1.7 


>0.8 


27.6±0.7 


27.8±0.7 


CL1358-z2<^ 






— 


>28.0 


>28.1 


CL1358-z3<i 


>2.2 


>2.0 


>1.2 


27.2±0.4 


>28.3 


A2218-zl«'f 


2.0±0.7 


>1.8 


>1.4 


26.7±0.2 


27.6±0.4 


A2218-Z2'' 


0.9±0.7 


>0.9 


0.5±1.0 


27.4±0.5 


26.9±0.4 


A2219-Z1 (lower left)s 




1.6±0.2 


0.5±0.3 


26.0±0.2 


26.0±0.2 


A2219-Z1 (upper right)^ 




>3.0 


1.0±0.5 


26.4±0.3 


26.3±0.2 


A2390-zl'^ 


>1.5 


>1.6 


0.0±0.4 


26.8±0.3 


26.6±0.2 


A2390-z2'^ 


>1.4 


>1.3 


-0.4±0.7 


27.4±0.6 


26.7±0.3 


A2667-zlf'i 


>2.1 


>1.4 


0.5±0.4 


26.6±0.3 


27.5±0.6 


A2667-Z2' 


2.1±0.7 


>2.0 


0.9±0.3 


25.6±0.1 


25.9±0.1 




z ~ 


9 J-drojjouts 








A2219-jlj 




2.0±0.3'' 


1.2±0.5'' 


26.7±0.2 


26.0±0.2 


A2667-jlj 


0.8±0.4'^ 


>1.3'' 


0.4±0.4'^ 


27.8±1.0 


26.4±0.2 



''Following the Richard et al. (2008) procedure, we have corrected the magnitudes we mea- 
sure in a 0.3"-diameter and 0.6"-diameter aperture by 0.3 mag and 0.6 mag, respectively, 
for the WFPC2+ACS and NICMOS data (see Appendix A). 

''This colour corresponds to r625 — Jiio for CL1358 and Abell 2218, V555 — ./no for Abell 
2390, and Vooe - JlW for Abell 2667. 

'^This colour corresponds to ^775 — Jno for CL1358 and Abell 2218, R702 — Jiio for Abell 
2219, and Isu - Jno for Abell 2390 and Abell 2667. 

''These sources are not detected at high significance (i.e., > 3it) in either the Jno or J/ieo 
bands in our NICMOS reductions. One possible reason that our calculated significance 
levels may be different from Richard et al. (2008) is that we account for sensitivity variations 
across the N1C3 detector in our weight maps (the sensitivities vary by factors of ~2 from 
region to region). 

''These candidates lie close enough to bright sources, as to make unambiguous detection 

and robust measurement quite difficult. 

'The measured ./no — Jfieo color here is much bluer than those found for 2 ~ 7 galaxies 
in the field (Bouwens et al. 2008) - strongly suggesting this source is not a 2 ~ 7 galaxy. 

SA2219-zl appears to be a blend of two sources (see Figure 1). While the lower left 
source is clearly a low redshift source, the upper right source seems to be a plausible 2 > 7 
candidate. 

''The 2 — J colours we measure for these sources appear too blue to correspond to dropout 
galaxies at 2 > 7. 

'These candidates have colours that are consistent with those of galaxies at 2 ~ 7. How- 
ever, it is difficult to be sure given the limited depth of the optical data over Abell 2667 
(Table 1: §4). 

J These sources are detected at >2<t at optical wavelengths in our reductions, strongly 
suggesting that they do not correspond to star-forming galaxies at 2 > 7. Richard et al. 
(2008) also note this fact and concede that these two sources are not particularly compelling 
2 > 7 candidates. 

''The colours tabulated here are relative to the if-band, not the Jno-band, and hence are 
0.6/i - Hieo, 0.8/i - i/ieo, and zsso - J^ieo- 



